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1 Introduction 



Lattice QCD aims to understand the strong interaction of hadrons from the first 
principles of QCD for quarks and gluons with the aid of numerical simulations. The 
physical quantities calculated with the method are many, ranging from the spectrum 
of light hadrons to a variety of electroweak matrix elements. Practical limitations such 
as finite lattice volume and spacing, and the use of the quenched approximation of 
00 \ ignoring sea quarks are gradually being lifted due to development of computer power, 

particularly those of dedicated parallel computers. In this talk we review progress in 

■ lattice QCD, focusing on efforts to calculate weak matrix elements relevant for the 
determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix . 

We start the review with discussion of the latest results for the light hadron 
spectrum and quark masses in sect. ||. In addition to being a topic basic to all of lattice 
Q_i| QCD, the recent progress in the quenched and full QCD calculations has yielded an 

o ■ interesting result, relevant for the CKM determination as well, that dynamical sea 

quarks have a significant effect on the value of the strange quark mass. 

The CP-violation in K meson decays continues to be an active topic with recent 
experimental reports on e'/£|]l|, Qj. In sect. we discuss the present status of the 
lattice calculation of the K meson mixing parameter and the decay amplitudes 

■ including the e'/e for which there has been recent renewed attempts. 
Understanding physics of B mesons is very important for the determination of the 

CKM matrix in the Standard Model and to detect the physics beyond it. It will be 
even more so as detailed data from the B factories, BaBar|| and Belle Q, become 
available. In sect. [| we present the lattice effort for calculating the weak matrix 
elements of heavy-light mesons. After a brief discussion on the theoretical issues of 
heavy quark calculations on a lattice, we present recent results on the B meson decay 
constant fs, the mixing parameter Bb and the form factors of semi-leptonic decays. 
For the decay constant serious attempts exist already to detect sea quark effects in 
its value. The other, more complicated quantities, are still calculated in the quenched 
approximation. For technical details and recent developments of heavy quark physics 
on the lattice we refer to recent reviews at the lattice conferences ||. 

In sect. H] we sketch to what extent the lattice results help constrain the CKM 
matrix determination. We conclude with a brief summary in sect. 0. 
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2 Light hadron spectrum and quark masses 



2.1 Light hadron spectrum 

Calculation of the light hadron spectrum has been pursued since 1981 because an 
agreement of the calculated spectrum with experiment provides a fundamental con- 
firmation for the validity of QCD in the non-perturbative low-energy domain, and 
also gives a measure of the reliability of the calculational techniques. 

A recent precision result of the CP-PACS Collaboration [|7| for the light hadron 
spectra in the continuum limit within the quenched approximation is shown in Fig. |l|. 
The up and down quarks are assumed to be degenerate {m u = = mi) in this 
calculation, and the experimental p and neutral 7r meson masses are used to fix the 
scale (lattice spacing a) and the light quark mass (mi). For the strange quark mass, 
two choices are compared, one employing the K meson mass (filled symbols ; K-input) 
and other with the <fi meson mass (open symbols; 0-input). Experimental values are 
shown by horizontal lines. 

This result shows an overall agreement of the light hadron spectrum in the quenched 
lattice QCD at a 5-10% level, as previously demonstrated by the GF11 Collaboration || 
However, it is also clear that there is systematic deviation between the quenched 
spectrum and experiment beyond the calculational error of 2-3%. In particular, the 
hyperfine splitting between the 0, K* meson masses and the K meson mass is smaller 
than the experimental one. 




Figure 1: Light hadron spectra in Quenched QCD. 
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Figure 2: and K* meson masses from K input as a function of a in full QCD (filled 
symbols), together with the quenched results (open symbols). 



As is shown in Fig. however, K* and (f) meson masses from K input agree much 
better in the 2 flavor full QCD|9|] than in the quenched QCD with the experimental 
values after taking the continuum limit. The effect of dynamical sea quarks is really 
important for reproducing the correct spectra. 

2.2 Light quark masses 

The masses of quarks are fundamental parameters of the Standard Model. Because 
of quark confinement, they can only be determined indirectly from a comparison of 
experimental hadron masses and their theoretical prediction in terms of quark mass 
parameters. For these reasons, much effort in lattice calculations have been devoted 
for their extraction. In this section we discuss the determination of strange quark 
mass m s which, as we shall see below, has an important impact on the e'/e in K 

meson decays. 

A summary of results for m^ IS (2 GeV) in quenched QCD with the Wilson quark 
action[10|, compiled by Bhattacharya and Gupta in 1997[|ll[], is shown in Fig. |3| as 



a function of a together with linear continuum extrapolations as explained in the 
caption. 

From this analysis for results with the Wilson quark action and similar ones for 
the other types of quark actions, they concluded as their best estimate of MS mass at 
2 GeV in quenched QCD to be mf^(2GeV) = 1 10(20) (ll)MeV, where the first error 



is estimated from the difference in the results among Wilson, Clover fll2|| and KSfl3|j 



quark actions, and the second one is the uncertainty of the scale (1/a) determination. 
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a (GeV -1 ) 



Figure 3: Strange quark mass in quenched QCD with Wilson quark actions. Two 
continuum extrapolations are shown for each input: one is the fit using all but the 
CP-PACS data, the other is the fit using only the CP-PACS data. 



This value is already close to the lower bound from QCD sum rules [H , mf IS (2 GeV) 
> 90 - 100 MeV. 

Two problems were manifest in their analysis. First, as can be seen from the 
figure, the large spread of data at finite a makes reliable continuum extrapolation 
difficult. This difficulty can be overcome by high precision calculations. Indeed the 
CP-PACS data in the figure is precise enough for a reasonable linear extrapolation. 

A more difficult problem, made quite clear by the CP-PACS precision data, is that 
the value of strange quark mass from .fT-input disagrees with the one from 0-input, 
even in the continuum limit. This is one of the manifestation of the quenching error. 

The full QCD result with 2 flavors of dynamical quarks from the CP-PACS col- 
laboration this year|| is shown in Fig. |j. Here the strange quark mass m^f s (2GeV), 
defined through the axial-vector Ward-Takahashi identity, is plotted as a function of 
a. Filled symbols in the figure are full QCD results obtained with the Clover quark 
action and renormalization group (RG) improved gauge action, while open ones are 
the previous quenched results with the Wilson quark action and the ordinary plaque- 
tte gauge action^. Results from both i-T-input (circles) and 0-input (squares) are 
given as in the previous figure. Numerically, we obtain in the continuum limit that 
mf^ (2 GeV) = 89.3 (6.6) MeV from if-input and 93.7 (8.8) MeV from 0-input in 
the 2 flavor full QCD, while 115 (2) MeV and 143 (6) MeV are the values for the K- 
and 0- input, respectively, in quenched QCD. 

It is very encouraging to observe that strange quark masses from two different 
inputs agree with each other within errors in the continuum limit. This reflects the 
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Figure 4: Strange quark mass in full QCD with Clover quark action (filled symbols) 
as well as the one in quenched QCD with Wilson quark action (open symbols), from 
K-input (circles) and 0-input (squares). 



fact, discussed already, that the hyperfme splitting among stranged mesons agrees 
well with experiment in the 2 flavor full QCD. 

It is surprising, however, that the quenching errors on the strange quark mass are 
larger than expected: 20% for iT-input and 40% for (p input. Moreover the strange 
quark mass in full QCD is much smaller than the previous quenched one, and is very 
close to the lower bound from QCD sum rules []L4j] . 

We quote the averaged value from the two inputs, 

mf^(2GeV) = 91(13)MeV, (1) 

as the CP-PACS result for the strange quark mass at 2 GeV in 2 flavor full QCD, 
which is equivalent to m^ 5 (m c ) = 100(15) MeV at the scale of the charm quark mass. 

Let us comment on the systematic uncertainties in this result: (i) The 2 fla- 
vor QCD simulation still neglects the dynamical effect of the strange quark itself. 
Although this effect is expected to be smaller than that from zero to 2 dynamical 
flavors, the value of the strange quark mass may be further reduced, (ii) The light- 
est dynamical sea quark masses used in the simulation corresponds the 7r over p 
mass ratio of m^/rrip = 0.55. One has to extrapolate the result obtained at these 
heavier quark masses to the physical light quark mass, where m^/rrip = 0.18. Their 
may be systematic uncertainties associated with opening of decay channels such as 
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p — > 2tt which are difficult to assess for the current range of heavy sea quarks, (iii) 
For the renormalization factor to convert the bare quark mass on the lattice into the 
renormalized MS quark mass in the continuum, the perturbative value calculated to 
one-loop order is used. It is certainly desirable to use the renormalization factor cal- 
culated non-perturbativelyljnj. For the present analysis, effects of higher order terms 
neglected are partly estimated by comparison of the quark mass evaluated through 
the axial-vector Ward-Takahashi identity and those from the mass parameter in the 
action. The influence on mf IS (2 GeV) is found to be 2-3 MeV, which has already 
been included in the estimate of the central value and the error, (iv) There exists 
also an uncertainty in the continuum extrapolation. Variations of the strange quark 
mass by a change of the ansatz in the continuum extrapolation, from linear in a to, 
for example, linear plus quadratic in a or a ■ aj^(l/a), are found to be 3-6 MeV, 
which have already been included in the error estimation. 



3 CP-violation in K meson decay 
3.1 Kq-Kq mixing parameter Bk 

The indirect CP- violation in K meson decays is parameterized by e, which is related 
to the CKM matrix element as 

r,[(l-p)A + B\B K = Q (2) 

where p and rj are Wolfenstein parameters, and A, B and C are numerical constants. 
Here Bk is the renormalization group invariant form of Bk, the Kq-Kq mixing pa- 
rameter, which is defined as the K meson matrix element of the AS = 2 four quark 
operator: 

R (K \sj^{l - ^ 5 )d ■ S7 M (I - 75)^0) ( v 

K = • {*) 

-{K \s^(l - 75 )d|0>(0| <hv(i - l,)d\K ) 

Thus the lattice calculation of Bk provides an important constraint on the CKM 
matrix. 

In Fig. H] we summarize the quenched lattice results for 5^(NDR, 2 GeV), renor- 
malized at 2 GeV in the MS scheme with the naive dimensional regularization (NDR), 
as a function of lattice spacing. These results have been obtained with a variety of lat- 
tice quark actions. The Wilson quark action PH [H| |I7] breaks chiral symmetry explic- 



itly, while the improved Wilson quark action (Clover or SW quark action) [|T2|, [18], [T 
reduces the magnitude of the chiral symmetry breaking. On the other hand, the 
Kogut-Susskind (KS) quark actionjl^] retains £7(1) subgroup of chiral symmetry, 



which is sufficient to guarantee the correct chiral behavior of the Bk parameter. 
Finally the do main- wall (DW) quark action P0|, recently developed and employed in 
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Figure 5: -B^NDR, 2 GeV) in quenched QCD with KS, Wilson/Clover and Domain- 
wall quark actions, as a function of the lattice spacing a. Two solid lines represent the 
continuum extrapolation by simultaneously fitting Bk from two different operators 
of the KS quark, with the constraint that two agree in the continuum limit. 



lattice QCD calculations [^1], 22|, maintains both chiral and flavor symmetries at the 
expense of introducing the extra dimension. After taking the continuum limit all 
formulation should give the same answer. 

As seen from the figure, the result from the KS fermion approach of the JLQCD 
Collaboration [p| is the most extensive and statistically accurate. They employ two 
different lattice operators for Bk, and matching to the continuum operator is made 
with the one-loop renormalization factor. In order to take into account the 2-loop 
ambiguity, in addition to the a 2 scaling violation of the KS quark action, they fit the 
results of Bk from two operators simultaneously by the form 

B K (a) = B K + b a 2 + c a^l/a) 2 , (4) 

with a unique value of Bk in the continuum limit. In the figure the fit is represented 
by solid lines. By this procedure the JLQCD collaboration obtain 

5^(NDR, 2GeV) = 0.628(42), (5) 

which corresponds to Bk = 0.87(6). The error in the continuum limit is much larger 
than those of individual data at non-zero a, due to the two-loop ambiguity a 2 in 
eq. (f§) . More accurate results requires precision determination of the renormalization 
factor in some non-perturbative way. 
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So far the degenerate quark masses, m s = nid, have to be used to avoid a quenched 
pathology that diverges in the limit that m s ^ — > 0. An 4 to 8% increase 
of Bk is predicted by chiral perturbation theory (%PT) for non-degenerate, physical 
quark masses p4l. As far as the quenching error is concerned, the value of Bk in 3 



flavor full QCD is found to increase by 5 % p5[. However this should be considered as 



a preliminary estimate, since no chiral and continuum extrapolation have been made 
in the calculation, and the degenerate quark masses are still used. 



3.2 Direct CP- Violation 



The direct CP-violation in K meson decays parameterized by e' is important since a 
non-zero value of e' strongly supports that the complex phase of the CKM matrix is 
the source of CP-violation in the Standard Model. Recent experiments consistently 
suggest that e'/e is non-zero, and the world average of experimental values is now 
21.3(2.8) xl(r 4 [|g. 

In order to see whether this value is consistent with the standard model prediction 
or suggests new physics, one has to theoretically calculate e'/ e within the framework 



of the Standard Model. An approximate formula takes the form[26 



e'/e ~ Im(A t ) x 13 



HOMeV 
m s (2GeV) 



x 



2.5 



MS I 

340MeV ) 



(6) 



where B, 



1/2(3/2) 
6(8) 



is the matrix element of the QCD (electroweak) penguin operator in 



the AI = 1/2(3/2) K meson decay, normalized by the vacuum insertion estimate, 
At = V t dV t * s , and Q v +r)' is the suppression factor by the isospin breaking effect in the 
quark masses. Since this formula indicates that the e' / e increases if the strange quark 
mass decreases, it is interesting to calculate its value using the lattice result of m s 
obtained in full QCD, which is smaller than the previous estimate. 

We have evaluated e'/e employing a more precise formula^] and taking the stan- 
dard estimate of parameters that B l J 2 = 1.0(3), B 3 8 /2 = 0.8(2), A§L = 340(50) MeV, 
ImAt = 1.33(14) x 10~ 4 and Q v +r]' = 0.25(8). For comparison we use two values of the 
strange quark mass, mf IS (m c ) = 100 (20) MeV from the lattice result and 130(25) 
MeV often employed in phenomenology. Here we conservatively add 5 MeV to the 
error of lattice result. 

We list the result in Table [I] varying the strange quark mass within the estimated 
error band; the other parameters are set to be their central values. Taking ambiguities 
in other parameters into account, the estimation in the Standard Model for e'/e, for 
a small strange quark mass indicated by the recent 2 flavor lattice QCD calculations, 
is still consistent with the experimental value. 
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m s (m c ) = 100(20) MeV: this talk 


m s (m c ) = 130(25) MeV:Ref. 26] 


m s (m c ) 


e'/e 


m s (m c ) 


e'/e 


80 MeV 
100 MeV 
120 MeV 


21.2 x 10" 4 
13.0 x 10" 4 
8.5 x 10" 4 


105 MeV 
130 MeV 
155 MeV 


11.5 x 10" 4 
7.0 x 10~ 4 
4.4 x 10" 4 



Table 1: The strange quark mass and e'/e 



It should be noted that only the strange quark mass is varied here while other 
parameters are kept fixed. This has inherent uncertainties since other parameters 
such as B^ 2 and B^ 2 also depend on the strange quark mass. (See ref. []27| for a 



quark mass independent parameterization of these matrix elements.) For a complete 

1 /2 3 /2 

prediction for e'/e from the Standard Model, we need to evaluate B G and B 8 also 
by lattice calculations. This possibility will be considered next. 

The direct calculation of the K — > tttt amplitude, in particular for the AI =1/2 
case, is notoriously difficult, mainly due to the euclidean nature of space-time on the 
lattice |28[]. Instead of the direct calculation, one may evaluate K — > tt amplitudes, 



which are related to the K — > tttt amplitude at the lowest order of x?T as follows p9| 



(tt + tt-\0,\K») = (vr+IO, - a i O sub \K + ) ■ ° K 2 ""\ (7) 

m MJM 



1 1 

m K — ml 
7T 



where Oi is a AS* = 1 four quark operator, and aij is determined by 

= (0\Oi - aiO sub \K°) (8) 

with O su b = (rrid + m s )sd + (m^ — m s )'s r y^d ) which can mix with Oi under operator 
renormalizations. Here tuk and m n are physical K and n meson masses, while rriM is 
an unphysical degenerate K and n meson mass, used in the lattice calculation of the 
K — > 7r matrix element. Therefore the left-hand side of eq. (0) should be independent 
on irtM at the lowest order of x?T. For this method to work, chiral symmetry of the 
lattice quark action is crucially important. 

Recently, applying this method with the KS quark action or the domain- 
wall quark action |p2[ |3jJ , both of which have good chiral properties, statistically 
meaningful results have been obtained for K — > titt amplitudes, relevant for the 
A = 1/2 rule. In Fig. || Re (A2)/Ke (Aq), where Aj is the K — > tttt amplitude for the 
final 7T-7T state with the isospin /, is plotted as a function of m 2 M (denoted as m\ in 
the figure), in both quenched and full QCD with the KS quark action[^D[. Although 
the ratio is still larger than the experimental value (AI = 1/2 enhancement is still 
smaller than the experimental one) and the strong m 2 M dependence, mainly caused 
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Figure 6: Re(A 2 ) /Re(A ) in quenched QCD with KS quark action at (3 — 6.0(crosses) 
and (3 = 6.2 (diamonds), as a function of the (unphysical) K meson mass squared, 
together with the one in full QCD at f3 = 5.7( suqares). A horizontal dotted line rep- 
resents the experimental value of the ratio, while the vertical one is the experimental 
value of K meson mass squared. 



by Re (A2), is observed, the result suggests that the method is useful. Calculations 
with better control over systematic errors as well as statistical errors will be the next 
step. A similar result for Re A with small m 2 M dependence is recently reported, by 
employing domain- wall quarks in the quenched approximation [F2^, [Hj]. 

Although results for e'/e have also been obtained by the same method with the KS 
quark action, it has been found that effects associated with 1-loop renormalization of 
weak operators are too large to obtain a reliable estimate of eyeflSOf. Partially employ- 



ing non-perturbative renormalization factors, a negative value of e'/e is obtained ||30||. 
Recently a large negative value of e'/e is also reported from the calculation by domain- 
wall quarks with non-perturbative renormalization factors p2|, |31j| . The value of B 6 



is also negatively large. However, further confirmations of these results, which in- 
clude a systematic study of scaling violations for the continuum extrapolation and 
the validity of the reduction method via x?T, will be needed for establishing a lattice 
estimate of e'/e. 

For the AI = 3/2 decay, on the other hand, the estimation by the lowest order 

3 /2 

of the ;\RT gives reasonable result. For example, the quenched results of B 7 and 
B^ 2 are given in Table |2|. Furthermore the direct calculation of the AI = 3/2 
K — > 7T7T decay amplitude successfully gives (n + n°\Q + \K + ) = (8.9-10.6) xl0~ 3 GeV 3 
in the quenched approximation [Q, in good agreement with the experimental value 
10.4xl0~ 3 GeV 3 . 
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References 


action 


a 


Renormalization 


B 3 7 /2 (2 GeV) 


B 3 8 /2 (2 GeV) 


GBS96[17| 

KGS97|32| 

Rome97|19| 

Rome97|19] 

UKQCD98|18| 


Wilson 
KS 

Clover 
Clover 
Clover 


-> 

7^0 
7^0 


Perturbative 
Perturbative 
Non- Pert urbat i ve 
Perturbative 
Perturbative 


0.58(2)(5) 

0.62(3)(6) 

0.72(5) 

0.58(2) 

0-58(i)(l) 


0.81(3)® 
0.77(4) (4) 
1.03(3) 
0.83(2) 
0.80(8)0) 



Table 2: Lattice estimates for 5 7 and B s in the quenched approximation. 



4 Weak matrix elements of heavy- light mesons 
4.1 Lattice action for heavy quarks 

A practical problem for lattice study of heavy quark physics is that the cut-off 1/a 
in the current simulations, which is typically 1 ~ 4 GeV, is smaller than the b quark 
mass, m\j ~ 4 GeV. Therefore the discretization error of order m^a is expected to 
be large. If the Wilson/Clover action is used for the heavy quark, the numerical 
simulations has to be restricted to the charm quark mass region, to avoid the large 
discretization errors, and the results has to be extrapolated to the b quark mass, 
with the help of results in the static limit. We call this method as the extrapolation 
method below. 

It has been pointed out ||34|| , however, that direct simulations at the b quark mass 
is possible if the large discretization errors is removed by a careful interpretation of 
the results. We call this the direct method. 

Although the continuum limit can be taken in both methods, the a dependence 
of the results is complicated and non-linear in the latter. 

If the non-relativistic formulation of quark is used instead, one can directly simu- 
late the b quark more easily and accurately. This method, called lattice NRQCD||35||, 



is an effective theory defined by an inverse expansion in the heavy quark mass. Tak- 
ing the continuum limit is not possible with this method, and the systematic error 
increases for the charm quark. 

Since there is no best method so far, one always has to check that different methods 
gives consistent results within estimated errors. 

4.2 Leptonic decay constants 

One of the simplest quantities in weak matrix elements is the leptonic decay constant 
of a heavy-light meson, which is important to determine V t d and V ts from neutral B 
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Figure 7: Jb, from different methods in quenched (open) and full (filled) QCD. 



meson mass differences: 

C /If 2 

AM q = -L^Lr lBq S{m t /M w )f B B Bq \V tq \\ (9) 

where denotes the decay constant of the B q meson. In particular a significant 
effort among lattice groups has been devoted to the calculation of fs = fB d , which 
has not been experimentally measured yet. 

After many years of calculations, results for fs from several groups with different 
methods are gradually converging. In Fig. |^, the summary of the latest results for 
f B form different methods, in both quenched (Nf = 0) and full (Nf = 2) QCD, are 
summarized. 

Results in quenched QCD, represented by open symbols, from the extrapolation 
method|§, NRQCD |37|, |38], g| and the direct method^ 0, @, are consistent 
within 10 % errors, showing that the lattice methods currently employed indeed work 
for b quark. Recent calculations in full QCD with 2 flavors of dynamical quarks (Nf = 
2), denoted by filled symbols, from NRQCD |43|, M and the direct method|(| g|, g5[, 
indicate that fs increases by 10-20 %. 

A summary of leptonic decay constants is given in Table || The present best 
estimate of fs in 2 flavor QCD is 210(30) MeV, which is larger than the previous 
quenched result. The ratio, fsJIsi which has less systematic uncertainty, on the 
other hand, remains unchanged from quenched QCD to full QCD. This quantity will 
become important if both AM^ and AM S are experimentally measured. 
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Quantity 


N f = 


N f = 2 


Experiments 


fa 


170(20) MeV 


210(30) MeV 




fs s 


195(20) MeV 


245(30) MeV 




fo 


200(20) MeV 




< 310 MeV 


fD s 


220(20) MeV 




241(21)(30) MeV 


IbJ Jb 


1.15(4) 


1.16(4) 




IdJ Id 


1.10(6) 







Table 3: Lattice estimates for leptonic decay constants in Nf = and 2. 



4.3 B meson mixing parameter B 



D 

As is seen in eq. @, the mixing parameters of neutral B q mesons B Bq (q = d,s), 
defined by 

B B M = 8 " 75)9 ' " lb)qW > (10) 

-(B q \bj,{l - T5 )g|0)(0|fo T ^l - l,)q\B Bq ) 

are important to determine Vtd and Vu from AM q , together with decay constants f Bq - 
Lattice results so far are obtained within the quenched approximation at fixed lat- 
tice spacings. A compilation of lattice results for Q Bb (aO = (&s(M p ) /a s (Ms)) 2///3 ° Bb(/j) 
at fi = 5 GeV is shown in Fig. |8], as a function of a heavy meson mass inverse, 
1/Mp. The coupling factor is introduced to cancel the ln(Mp) dependence in the B 
parameter . 

Data in the figure are classified into 3 distinct groups. The first is the calculation 
from the static-light approximation, where the heavy meson mass is fixed to be infi- 



nite. The results are consistent among the Kentucky group[47], Gimenez- Reyes [48. 49 
and UKQCD0], though the error of the Kentucky group is rather large due to the 
ambiguity of 1-loop renormalization in Wilson quark action. Therefore, we take an 
average over results of the latter two groups, who employ the Clover quark action, as 
the current lattice estimate, giving B Bd (mb) = 0.80(5) (B Bd = 1.28(8)). 



The second result comes from a NRQCD-light method by the Hiroshima group [ 52 1 . 
They observe a heavy meson mass dependence: B B decreases as B meson mass de- 
creases from the static limit. As is seen from the figure, their result in the static limit 
is reasonably consistent with the previous result from static-light calculations. 

There are several calculations using the extrapolation method. They are rather 



old results from Wilson quark action|53], [54|, |17 , except UKQCD's this year with 



the Clover action[Q. Although the results are internally consistent, they seem to 
differ from the static result in the static limit and the NRQCD result at the B meson 
mass, as is seen in the figure. More work is needed to reduce this uncertainty and to 
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Figure 8: A compilation of quenched results for $_b s (5 GeV) as a function of a heavy 
meson mass inverse, 1/Mp. 



obtain the reliable estimate for B B . 

We take the static result as our tentative estimate for B B , since systematic errors 
are best controlled in this method so far: 



B Bd (m b ) = 0.80(5), (B Bd = 1.28(8)), £ ad = hs ^ s = 1.17(6). ( 1 1 ) 



fB d \]B Bd 

Combining this value with the previous estimate of f B , one obtains 



4th d f Bd = 240(36)MeV. 



(12) 



4.4 Form factors of semi-leptonic decays 

In this section form factors of semi-leptonic decays of B or D mesons are considered. 
For pseudo-scalar to pseudo-scalar decays there are two form factors, f + (q 2 ) and 
f°(q 2 ), which are defined as 



(P(k)\V»\H(p)) = f + (q 2 



(p + ky 



m P u 



(13) 



here q = p — k is the momentum transfer, H = B or D, and P = D, K, 77 or ir. Four 
form factors for pseudo-scalar to vector decays, V(q 2 ), A 12 (g 2 ) and A(q 2 ) are defined 



as 



(V(k,e)\V»\H(p)) = -ffl-e^^ej 



14 



(V(jfe,e)|A1#(p)) = i(m H + m v )A 1 (q 2 ) e*" -i A ^ e* ■ p (p + kf 



m H + my 



+t^r 1 2m v e* -p 



where V = D*, K*, <\> or p. 
4.4.1 D -»■ ifW^i/, 7r(p)£z/ 

We first consider D meson decays to light mesons, K, K*, ir or p. The purpose of this 
calculation is two- fold: Assuming that V cs = 1 — A 2 /2, one can establish the validity 
of lattice calculations for semi-leptonic form factors. Conversely one may extract V cs 
or V c b directly from experimental data, using the lattice estimate of the form factors. 

The calculation of these form factors is easier than others, because direct simula- 
tions at m c or m s can avoid subtlety of the quark mass extrapolation. In addition the 
fact that all physical as well as unphysical q 2 regions are covered makes interpolation 
to q 2 = possible and easy. 

As an example, in Fig. |9| a lattice result [57|], already in 1995, for the form factor 



f + (q 2 ) of the D —>■ Klv decay, is given as a function of q 2 in quenched QCD. As this 
figure shows, lattice methods work well for these decays. For a summary of quenched 
estimates for these form factors at q 2 = 0, see Ref.||58||. 



1.5 
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0.5 



-0.1 0.0 0.1 



0.2 



Figure 9: The form factor f + (q 2 ) for D — > Klv decay, as a function of q 2 
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Figure 10: \J r B _ tD *{uj)V c b\ as a function of u from DELPHI collaboration. 



4.4.2 B -> 

The large branching ratio of B — > D^£u decays makes a determination of |V^| very 
precise. Differential decay rates are related to |K&| 2 and the form factor Tb^d by 

— {B -> I)W&/) = (known factor) x |V r J 2 |.F B _> £ K.)M| ( 14 ) 
do; 

where u = Vb • v D (*) is the velocity transfer. For example, the very accurate result 
for \J 7 B-*D*{w)V c b\ from the DELPHI collaboration is shown in Fig. [HJ Taking the 
u -> 1 limit, one obtain I-Fb^m (1)1^1 = (37.95 ± 1.34 ± 1.59) x 1(T 3 [|1|. 

Now the most important task of lattice QCD is the determination of |JF B ^ I) (*)(1)|, 
in order to extract |V^| from experimental results. They are given by 



Fb^b(1) = Ml)- W \ mJ L(l) (15) 

?™B + ?™D 

Fb^b*(1) = Mi), ( 16 ) 

in terms of new form factors, defined by 

= M«)<« + «/), + *-(«)<«-«%] (it) 

(D*(v')\^ l5 b\B(v)) 



y/m B m* D 



(u + l)elh Al (u). (18) 



Since the heavy quark symmetry implies JFb^d(I) = Fb^d*{\) = 1 in the m&, 
m c — > oo limit, the heavy quark mass dependence has to be determined precisely. 

These form factors have been evaluated recently by the Fermilab group |)9|, using 
the direct method with Clover quark action, in quenched approximation at a -1 ~ 
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Lattice 


QCD sum rule 




1.058 (16) (+^ 4 ) 
0.935 (22) (±|I) 


0.98(7) 

0.91(6)/0.92(4) 



Table 4: Lattice results for \J 7 b-*d*(^)\, together with those from QCD sum rule. 



1 GeV. In order to cancel systematic as well as statistical errors of lattice matrix 
elements, they employ the ratio method, through which |/i + (l)| 2 , for example, is 
obtained as 

\ h B^p m] 2_ (D\VS b \B)(B\V^\D) f . 

1 + {)l (D\V<r\D)(B\V bb \By K ' 

The precison of the result makes it possible to determine the heavy quark mass 
dependence of h's up to 1/rnf c for h + (l) and l/ml c for ^^(l). 

The summary of their results, together with results from the QCD sum rule 
approach ||60||, is given in Table |j. Already errors of lattice results are comparable 



to or even smaller than those of QCD sum rule results. Since the deviation from 1, 
not the value itself, is important for F B ^, D {*)(1), the systematic as well as the sta- 
tistical errors should be reduced further. Improved calculations with smaller lattice 
spacings for the continuum extrapolation and/or the inclusion of the full QCD effect 
should be performed as the next step. Results from such calculations, together with 
experimental data, will reduce errors of \V c b\ significantly. 

As far as the shape of F B _^ D (») (u) is concerned, the lattice calculation has a longer 
history than that for F B _ t , D {*)(l). It is important to calculate the shape of form factors, 
as a check of theoretical method and for a reduction of errors in the extrapolation of 
experimental data. 

Heavy quark symmetry implies that these form factors are related to one universal, 
mass-independent function, £(u;), the famous Isgur-Wise function: 

h+, Al (u) = (1+ 7 W^)+/3+,aH)£H, (20) 

h-(u) = (t-M + zM^H, (2i) 

where /3± tAl represent radiative corrections, while j± jAi (uj) represent heavy quark 
mass dependences. The Isgur-Wise function has already been calculated on 
the lattice. For example, the recent result for h+(u)/(l + (3 + (u)) from the UKQCD 
collaboration |nj with the Clover quark action is given in Fig. [TTJ. They employ four 



different heavy quark masses, and find the mass-independence of h + (uS)/(l + (3 + (uj)), 
which suggests that the correction to the heavy quark mass limit is indeed small. 
Therefore, within errors, h + (u)/(l + (5 + {uj)) can be identified with the Isgur-Wise 
function, £(u). In future investigations it will be important to extract mass dependent 
term, 7± J i 1 (w), in order to construct the shape of F b _+ D (*)(uj) precisely. 
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Figure 11: The radiatively corrected form factor h + (u) in quenched QCD from 
UKQCD collaboration. Different symbols correspond to different heavy quark masses, 
while the light quark mass is kept fixed. 



4.4.3 B -> nlv, piv 

The last quantity of our concern is the form factors of B — > ir, p decays, which are 
rare decays, and can be used to determine \V u i\. 

Lattice calculations for these form factors are restricted to large q 2 regions only, 
since momenta of B and n/p have to be smaller than 1/ct, to avoid large discretization 
errors. Therefore experimental data for the partial decay rate in the large q 2 region 
are necessary for such lattice calculations being useful. Indeed statistics in the recent 



experiment is precise enough to attempt an initial study for the partial decay rate|62 
and the quality of data is expected to be further improved. 

The UKQCD collaboration calculates the form factor relevant for the differential 
decay rate of B — > piv in the quenched approximation at a ^ 0, using the Clover 



quark action |63|. As shown in Fig. [12], fitting lattice data near q = g^ ax by the form, 
dT(B piv) 10 12 9/ 9 9 .. . , 

— —T% ITTl^ = c 1 + b y " W x Phase space factor . 22 

dq 2 \V ub \ 2 

they obtain 

c = 4.6 ± 0.7GeV b = (-8^) x 10~ 2 GeV" 2 . (23) 

As the precision of the differential decay rate will be further improved in future 
experiments, this information will be useful in the extraction of \V u b\. 
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Figure 12: The form factor B — > p£t/ decay as a function of g 2 in quenched QCD 
from the UKQCD collaboration. The solid line represents the fit given in the text, 
and vertical dotted line shows q^ ax for B —> D* decay. 



There are two form factors for B — > islv decay, defined as 

2 2 1 11 

Tn ~ — Tn „ TP ^ 

( 7 r(A ; )i^|5(p)) = / + (g 2 



2 2 

777.7, — ml 



+ W) ™b™^. (24) 



g 2 



2 (fi^7r^)cx|K b ri/ + (? 2 )| 2 , (25) 



Only f + (q 2 ) contributes to the differential decay rate as 

dr 

dq 

while the contribution of f°(q 2 ) is negligible in the decay rate since q^L^ oc mi, where 
L M is the leptonic weak current and mi is the lepton mass. However f°(q 2 ) may be 
used to check lattice calculations by the relation that /°(q , m ax ) = }'b/ fn, implied by 
the soft pion theorem in the m n — > limit. 

It has been found that the calculation of these form factors is rather difficult. 
The chiral extrapolation of f°{q 2 ) linear in the light quark mass does not appear to 
satisfy the relation of the soft pion theorem^, Although it is claimed that the 
inclusion of the a square root term in the chiral extrapolation solves the problem |65|, 



at this moment we have to conclude that lattice results are premature for a detailed 
comparison with experimental data, and that more work will be needed on this form 
factor. 



5 Impact on determination of CKM matrix 

Let us collect the lattice results and discuss their impact on the CKM matrix deter- 
mination. We employ the four standard constraint relations given below to determine 
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Experimental inputs 


Theoretical inputs 


quantities 


value 


quantities 


Buras99 


Lattice 


\Vcb\ 


0.040(2) 


B K 


0.80(15) 


0.87(6) 


\v ub \ 


3.56(56) • 10~ 3 


i B 




210(30) MeV 


\e\ 


2.280(13) • 10~ 3 


B B d 




1.28(0.08) 


AM, 


0.471(16) ps- 1 


\jB~BfB 


200(40) MeV 


240(23) MeV 


AM, 


> 12.4 ps- 1 


t 


1.14(8) 


1.17(6) 


m t 


165(5) GeV 









Table 5: Input parameters. 



the unitarity triangle (see ref. 



p 2 + rj 2 



\ 10 a 



26J for detailed notations). 

A 2 , 



1 



V, 



ub 



v, 



cb 



1 — p) 2 + rj 2 



[1 - 
V td 



p)A 2 rj 2 S (x t )+P (e: 



A 2 B 



K, 



with AM d 



"7, oVb ■ m B 

D7T 



B B jB d m 2 w S (x t ) \V t 



td\ 



l-p) 2 + 7] 2 > ^ 



10.2/ps 



\ (AM, 



(26) 
(27) 

(28) 
(29) 



where p = p(l — A 2 /2) and r\ = r](l — A 2 /2). Experimental inputs, which come from 
b — » u transition for eq. (|26|) , K — K (e) for eq. (p7|) , and B — B mixings for 
are summarized in Table El. As theoretical inputs, we consider two 



eqs. (|2|, |2|) 

choices, also specified in Table [5|: one choice based on the lattice results summarized 
in this review, and the other choice from Ref. |[26|| . We shall refer to the latter as a 
standard one. 

In Fig. [13| the plot on the left shows the output in the p— rj plane from the standard 
input. The shaded region is allowed within 1 a of each constraint. We take a flat 
distribution of each error in this analysis. In the plot on the right the output from the 
lattice input is represented by thin curves (denoted as "present"). Due to the larger 

value of \j B B f B than the one in the standard input, more positive values of p are 
favored. Thick curves (denoted as "ideal" ) are constraints from lattice inputs without 
errors. In both plots the constraint from AM, is represented by dashed curves. 

From these figures we conclude the following. For a better determination of CKM 

matrix, lattice calculations have to reduce errors in y B B f B and precisely determine 
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the central value of both and y Bgfs, with less systematic uncertainties (in the 
continuum limit of full QCD). On the other hand, experiments need to reduce errors 
in \V c b\ and \V u b\, where lattice calculations can also contribute, and measure AM S 
instead of a lower bound available at present. 

6 Summary 

In this review we have presented the status of lattice QCD calculations of the quan- 
tities relevant for further understanding of the Standard Model and its limitations. 
The main results may be summarized as follows. 

• In the continuum limit of full QCD with Nf = 2, lattice QCD predicts that 
m^ IS (m c ) = 100 (15) MeV, which is smaller than expected. This small value of 
the strange quark mass tends to increase the estimate for e'/e in the standard 
analysis. 

• The continuum limit of quenched calculations for the Kq — Kq mixing parameter 
gives Bk = 0.87(6). Clearly the full QCD estimate is called for. 

• Calculations in full QCD with Nf = 2 for the B meson decay constant estimate 
fs = 210(30) MeV , which is found to be larger than previous quenched results. 

• The Bq — Bq mixing parameter I?B(m&) = 0.80(5) in the static limit of the 
quenched approximation. The next step is to estimate the heavy meson mass 
dependence. 
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• Lattice methods for form factors of D — > K^*\ n, p decays are now well estab- 
lished. There is a promising method for the calculation of J r B _ fD (*)(l) to extract 
\V c b\- More works are still needed for B — > 7r (p) decays. 

As a conclusion we stress that systematic investigations for all these quantities in the 
continuum limit of full QCD will be the next target. 

I thank Drs. P. Collins, S. Hashimoto, H. Shanahan, and A. Ukawa for informative 
communications and useful discussions. 
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Discussion 



Masanori Yamauchi (KEK): Do the errors quoted represent one standard devi- 
ation, just like experimental statistical errors? 

Aoki: Yes, the quoted numbers generally represent one standard deviation of sta- 
tistical errors, while some of them also include systematic errors. 

Jeff Richman (UCSB): I would like to comment that there is no fundamental 
obstacle to measuring the rate for B — > piv at high q 2 . In fact, the experimental 
acceptance is good in this region, and CLEO has now obtained a crude measurement 
of dT I dq 2 for B — > p£v, which is limited only by statistics. If we meet again in three 
years, perhaps the lattice will have a solid prediction for the rate at high q 2 and 
experimenters will have a solid measurement. 

Aoki: I agree with you. I think that lattice theorists and experimentalists in this 
field should communicate with each other. 

George Hou (National Taiwan University): The strange quark mass m s not 
only enters e'/e, but it also enters rare B decays. In the CLEO fit discussed in 
Poling's talk that obtained a somewhat different value for 7 (or 2 ) than usual, one 
also obtains a m s value consistent with the lower value that you have reported. 

Aoki: That is an interesting information. I would like to see the detailed analysis 
in the future. 

F. Constantini (Pisa University): It is true that the decrease of the strange q 
mass from 130 MeV to 100 improves the agreement between the lattice result and the 
experimental value of e'/e. However, to agree with the present experimental values 
you need to go as low as 70 MeV. Moreover the dependence on the strange quark 
mass acts on the numerator of the expression and not only on the denominator. 

Aoki: Since there exist ambiguities also in other parameters such as B$, B$ and 
A.jjg, we do not need the strange quark mass as low as 70 MeV. The explicit strange 
quark mass dependence in the numerator of the expression is included in my estimate. 
As mentioned in the talk, however, the implicit strange quark mass dependence in 
B 6 and B$ is not considered here. 

Matthias Neubert (SLAC): In fact the quantity e'/e does not depend at all on 
the strange quark mass. It enters only by convention when you define B 6 . It seems 
to me you should not use a lattice estimate of m s unless at the same time you use a 
lattice estimate of B 6 . 
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Aoki: Rigorously speaking, I agree with you. According to this criteria, however, I 
have to say that there is no reliable theoretical estimate for e'/e. Now it seems for 
me that only lattice calculations can give a reliable theoretical estimate for e'/e, and, 
as I mentioned in my talk, serious efforts on such calculations have just started. 

B. F. L. Ward (University of Tennessee): Could you comment on the prospects 
for calculating the rate of B — > tttt on the lattice? 

Aoki: There are two main difficulties for calculating the rate of B — > tttt on the 
lattice. One is the problem of dealing with decays into 2 particles suck as K — > tttt 
decays, the other is the difficulty of putting a heavy quark on the lattice. Therefore, 
no lattice calculation has been made for the rate of B — > tttt so far. 

Helen Quinn (SLAC): When you put your errors on lattice quantities for which 
you have only a quenched calculation I believe that these errors include statistics and 
errors due to extrapolation but no estimate of the errors due to quenching. Is this 
correct, and if so how can one best estimate the additional uncertainties due to the 
quenched approximation? 

Aoki: Yes. In general the errors due to quenching are not included, and it is 
difficult to estimate them unless the corresponding full QCD calculation is really 
performed. One might naively think that the quenching errors are 10-20%. In some 
calculations, the quenching errors are partly estimated by the variation of the result 
under the change of experimental input parameters. For example, the difference of 
the strange quark mass between X-input and $-input, which is about 20-25%, may 
be interpreted as the quenching errors. However this kind of the estimation is far 
from being quantitative. 



26 



